The zebrafish transcriptome during early development

<p>Abstract</p> <p>Background</p> <p>The transition from fertilized egg to embryo is accompanied by a multitude of changes in gene expression, and the transcriptional events that underlie these processes have not yet been fully characterized. In this study RNA-Seq is used to compare the transcription profiles of four early developmental stages in zebrafish (<it>Danio rerio</it>) on a global scale.</p> <p>Results</p> <p>An average of 79 M total reads were detected from the different stages. Out of the total number of reads 65% - 73% reads were successfully mapped and 36% - 44% out of those were uniquely mapped. The total number of detected unique gene transcripts was 11187, of which 10096 were present at 1-cell stage. The largest number of common transcripts was observed between 1-cell stage and 16-cell stage. An enrichment of gene transcripts with molecular functions of DNA binding, protein folding and processing as well as metal ion binding was observed with progression of development. The sequence data (accession number ERP000635) is available at the European Nucleotide Archive.</p> <p>Conclusion</p> <p>Clustering of expression profiles shows that a majority of the detected gene transcripts are present at steady levels, and thus a minority of the gene transcripts clusters as increasing or decreasing in expression over the four investigated developmental stages. The three earliest developmental stages were similar when comparing highly expressed genes, whereas the 50% epiboly stage differed from the other three stages in the identity of highly expressed genes, number of uniquely expressed genes and enrichment of GO molecular functions. Taken together, these observations indicate a major transition in gene regulation and transcriptional activity taking place between the 512-cell and 50% epiboly stages, in accordance with previous studies.</p

Identification of Novel Transcribed Regions in Zebrafish (Danio rerio) Using RNA-Sequencing

Zebrafish (Danio rerio) has emerged as a model organism to investigate vertebrate development and human genetic diseases. However, the zebrafish genome annotation is still ongoing and incomplete, and there are still new gene transcripts to be found. With the introduction of massive parallel sequencing, whole transcriptome studies became possible. In the present study, we aimed to discover novel transcribed regions (NTRs) using developmental transcriptome data from RNA sequencing. In order to achieve this, we developed an in-house bioinformatics pipeline for NTR discovery. Using the pipeline, we detected 152 putative NTRs that at the time of discovery were not annotated in Ensembl and NCBI gene database. Four randomly selected NTRs were successfully validated using RT-PCR, and expression profiles of 10 randomly selected NTRs were evaluated using qRT-PCR. The identification of these 152 NTRs provide new information for zebrafish genome annotation as well as new candidates for studies of zebrafish gene function.Peer reviewe

The Zebrafish Orthologue of the Dyslexia Candidate Gene DYX1C1 Is Essential for Cilia Growth and Function

Peer reviewe

Differences in Gene Expression between Mouse and Human for Dynamically Regulated Genes in Early Embryo

Peer reviewe

Globin mRNA reduction for whole-blood transcriptome sequencing

The transcriptome analysis of whole-blood RNA by sequencing holds promise for the identification and tracking of biomarkers; however, the high globin mRNA (gmRNA) content of erythrocytes hampers whole-blood and buffy coat analyses. We introduce a novel gmRNA locking assay (GlobinLock, GL) as a robust and simple gmRNA reduction tool to preserve RNA quality, save time and cost. GL consists of a pair of gmRNA-specific oligonucleotides in RNA initial denaturation buffer that is effective immediately after RNA denaturation and adds only ten minutes of incubation to the whole cDNA synthesis procedure when compared to non-blood RNA analysis. We show that GL is fully effective not only for human samples but also for mouse and rat, and so far incompletely studied cow, dog and zebrafish.Peer reviewe

Novel PRD-like homeodomain transcription factors and retrotransposon elements in early human development

Transcriptional program that drives human preimplantation development is largely unknown. Here, by using single-cell RNA sequencing of 348 oocytes, zygotes and single blastomeres from 2- to 3-day-old embryos, we provide a detailed analysis of the human preimplantation transcriptome. By quantifying transcript far 50-ends (TFEs), we include in our analysis transcripts that derive from alternative promoters. We show that 32 and 129 genes are transcribed during the transition from oocyte to four-cell stage and from four-to eight-cell stage, respectively. A number of identified transcripts originates from previously unannotated genes that include the PRD-like homeobox genes ARGFX, CPHX1, CPHX2, DPRX, DUXA, DUXB and LEUTX. Employing de novo promoter motif extraction on sequences surrounding TFEs, we identify significantly enriched gene regulatory motifs that often overlap with Alu elements. Our high-resolution analysis of the human transcriptome during preimplantation development may have important implications on future studies of human pluripotent stem cells and cell reprograming.Peer reviewe

Characterization and target genes of nine human PRD-like homeobox domain genes expressed exclusively in early embryos

PAIRED (PRD)-like homeobox genes belong to a class of predicted transcription factor genes. Several of these PRD-like homeobox genes have been predicted in silico from genomic sequence but until recently had no evidence of transcript expression. We found recently that nine PRD-like homeobox genes, ARGFX, CPHX1, CPHX2, DPRX, DUXA, DUXB, NOBOX, TPRX1 and TPRX2, were expressed in human preimplantation embryos. In the current study we characterized these PRD-like homeobox genes in depth and studied their functions as transcription factors. We cloned multiple transcript variants from human embryos and showed that the expression of these genes is specific to embryos and pluripotent stem cells. Overexpression of the genes in human embryonic stem cells confirmed their roles as transcription factors as either activators (CPHX1, CPHX2, ARGFX) or repressors (DPRX, DUXA, TPRX2) with distinct targets that could be explained by the amino acid sequence in homeodomain. Some PRD-like homeodomain transcription factors had high concordance of target genes and showed enrichment for both developmentally important gene sets and a 36 bp DNA recognition motif implicated in Embryo Genome Activation (EGA). Our data implicate a role for these previously uncharacterized PRD-like homeodomain proteins in the regulation of human embryo genome activation and preimplantation embryo development.Peer reviewe

Prediction and validation of 4 randomly selected NTRs.

<p>The RNA samples used for validation were extracted from 50% epiboly. RNA-seq: RNA-Seq tracks (red boxes) based on the pooled RNA-Seq data; NTRs: Predicted structures of putative NTRs by our pipeline. Blue boxes represent fragments; Amplicons: Validations of the predicted NTR structures by RT-PCR. “Chr” indicates chromosome. A. NTR50; B. NTR88; C. NTR103; and D. NTR145.</p

Systematic workflow for the identification of NTRs using RNA-Seq datasets.

<p>A. RNA-Seq reads form fragments (shown in green, blue, red and purple), which further generate clusters for NTR identification. D1: the distance between a putative NTR and any annotated transcribed regions; D2: the distance between two putative NTRs; B. Systematic workflow of NTR identification. The numbers for fragments are circled in blue, while numbers for clusters are in red circles. Singleton fragments are one-fragment clusters with length over 50 bp.</p